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ABSTRACT 


A complete  data  processing  package  for  navigation 
and  total  geomagnetic  field  anomaly  is  described  for 
the  Hewlett  Packard  9830  calculator.  These  programs 
were  designed,  written  and  used  during  a two-week  site 
survey  funded  through  the  International  Program  of  Ocean 
Drilling  (IPOD) . These  programs  permit  1)  entry  of 
navigation  and  total  field  geomagnetics  data;  2)  data 
reduction  of  total  field  to  magnetic  anomaly;  3)  data 
editing,  and  4)  data  display  in  the  form  of  annotated 
ship's  track  and  profiled  data  along  the  ship's  track. 


INTRODUCTION 


A two  week  site  survey  funded  through  the 
International  Program  of  Ocean  Drilling  (IPOD)  took 
place  from  1 to  15  September  1976  aboard  the 
U.S.N.S.  LYNCH  (T-AG0R7).  Two  alternate  sites 
isochronous  to  the  Late  Mesozoic  site  AT2 . 1 (formerly 
termed  "2A")  were  surveyed.  It  is  the  purpose  of 
this  report  to  describe  the  data  processing  software 
written  for  the  Hewlett-Packard  9830  calculator 
available  on  the  LYNCH. 

There  was  a need  to  chart  the  navigation,  depths  and  magnetics 
data  soon  after  collection,  in  order  to  make  decisions  on 
subsequent  ship's  tracks.  The  W.H.O.I.  data  acquisition 
system,  general  purpose  computer  and  processing  software 
were  not  available  for  this  cruise.  Although  the  underway 
data  were  digitized  by  hand  and  key-punched  daily 

aboard  ship  onto  computer  cards  for  later  processing  at 
W.H.O.I.,  we  needed  an  alternate  processing  package  at 
sea . 


These  programs  were  written  at  sea  during  the  first 
week  of  the  cruise.  They  were  written  in  BASIC  and 
operated  within  the  capabilities  of  the  HP  9830.  The 
software  provided  the  following  capabilities: 

1.  Navigation  data  logging,  and  calculation  of 
speed  and  heading  between  fixes. 

2.  Even  five  minute  geomagnetics  position  interpolation 

3.  Total  magnetic  regional  field  interpolation, 
measured  magnetic  field  logging  and  anomaly 
calculation . 

4.  Edit  routine  for  navigation  data. 

5.  Edit  routine  for  magnetic  data. 

6.  Mercator  charting  with  time  annotated. 

7.  Mercator  charting  with  magnetic  anomaly  plotted 
as  a profile  parallel  to  the  ship's  track. 

Program  write-ups  and  listings  appea _ 
of  this  report. 


in  Appendices  4 to  9 
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EVOLUTION  OF  DATA  PROCESSING  SCHEME 


At  the  start  of  the  cruise  the  navigation  data, 
consisting  mostly  of  satellite  fixes  and  dead  reckoned 
positions,  were  hand  plotted  onto  a large  scale 
(16  inches  per  degree)  predrawn  Mercator  base  map. 

The  depth  data  were  transferred  by  hand  to  a similar 
map.  Since  the  depths  were  digitized  at  about  every 
five  minutes,  a program  was  written  to  calculate  and 
list  even  five  minute  interpolated  positions  between  the 
fixes.  A separate  program  was  written  to  calculate  the 
speed  and  heading  between  fix  pairs  and  was  used  to 
verify  the  positions.  Another  program  was  written  to 
accept  the  interpolated  five  minute  positions  and  the 
total  measured  geomagnetic  field  value.  It  calculated 
the  regional  field  at  each  of  the  five  minute  values  by 
linearly  interpolating  within  the  deqree  square,  given  the 
regional  geomagnetic  field  at  the  time  of  the  cruise  for  each 
corner  of  the  degree  square.  A listing  of  time,  position 
and  anomaly  were  then  made.  Each  step  of  the  processing 
required  a separate  entry  for  the  navigation  information. 
These  tasks  were  repetitious,  tedious  and  unnecessary. 

We  needed  a simple  but  unified  approach  to  the  data 
collecting  and  display. 

The  Hewlett  Packard  9830  calculator  on  the  ship  had 
limited  input/output  capabilities:  one  built-in  cassette 
read/write  unit.  Merging  of  data  between  separate  input 
devices  was  therefore  not  possible.  However,  the  task 
of  typing  over  5000  numbers  per  day  into  the 
calculator  was  an  incentive  for  improving  the  existing 
scheme . 

The  processing  system  that  evolved  uses  the  calculator 
memory  to  store  approximately  sixteen  hours  of  navigation 
and  magnetics  data.  The  memory  served  as  the  random 
access  storage  device.  Because  of  the  calculator  design, 
the  data  stored  in  COMMON  was  accessible  to  several 
programs  having  the  same  COMMON  specification.  The 
cassette  unit  provided  storage  for  the  data  already 
processed  and  allowed  for  program  overlaying. 
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Our  initial  goal  was  to  eliminate  the  redundant  task 
of  keying  in  the  navigation  data  three  times.  This  was 
accomplished  very  quickly.  Of  course,  once  the  data 
were  collected  and  properly  organized  in  memory,  many 
other  options  became  possible.  In  particular,  since 
the  calculator  was  connected  to  a flat  bed  plotter  we 
were  able  to  display  the  data  on  charts.  The  10"  x 15" 
plotter  was  quite  small  for  the  expanded  (16  "/degree) 
longitude  Mercator  scale  used  for  mapping  the  area, 
and  several  charts  were  required  to  cover  the  survey 
area.  These  'chartlets'  were  traced  onto  the  base  map. 

This  was  an  improvement  to  the  use  of  dividers  and  a 
Gerber  for  the  576  points  of  navigation  and  magnetics 
per  day. 

Three  programs  were  used  to  make  a "merged"  data 
file  consisting  of  data  value,  position  and  time  of 
observation.  The  first  program  accepted  the  fixes,  through 
the  keyboard,  and  saved  them  in  memory.  After  all  the 
fixes  were  keyed  in,  the  program  produced  a listing  of 
speeds  and  headings  between  fixes.  Missed  fixes,  mistaken 
data  entries  and  bad  dead  reckoned  positions  could  be 
identified  in  this  way.  An  edit  program  was  run  to  correct 
these  errors.  This  editing  routine  allowed  for  modifying, 
inserting  and  deleting  of  fixes.  After  errors  were 
removed,  the  third  program  was  run.  This  program  calculated 
even  five  minute  positions  between  fixes.  Two  hundred  of 
these  times  and  positions  could  be  stored  in  memory  at  a 
time.  These  data  could  then  be  stored  on  cassette  for 
later  use.  The  anomaly  calculating  program  used  this 
navigation  data  in  calculating  the  regional  field  values. 


PROGRAM  LIMITATIONS  AND  RESTRICTIONS 


There  are  several  limitations  in  the  programs  due 
largely  to  the  dispatch  with  which  they  were  w-itten. 

The  site  survey  took  place  in  the  area  bounded  by  30  N, 

20  N,  75  W and  65  W.  Data  entry  verification  was  much 
simplified  by  taking  advantage  of  this  fact.  Also,  the 
cruise  took  place  wholly  within  the  month  of  September. 
Consequently,  time  interval  calculations  could  be  much 
simplified.  The  programs  did  not  have  to  take  leap  years 
or  changes  in  month  into  account. 
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Care  was  taken  in  the  coding  for  error  checking 
at  data  entry.  Here  too  advantage  was  taken  of  the 
limited  geographic  scope  of  the  survey.  The  regional 
magnetic  field  values  based  on  the  1975.0  spherical 
harmonic  coefficients  have  a positive  gradient  in  the 
north-westerly  direction.  Whenever  four  regional  field 
values  were  entered  they  were  checked  for  this  trend. 
These  checks  of  the  operator  inputs  repeatedly  saved 
time . 

These  geographic  restrictions  will  have  to  be 
revised  in  order  to  use  this  processing  scheme  at 
other  survey  sites.  The  ease  of  editing  BASIC  routines 
on  the  HP  9830  calculator  makes  such  input  checks 
readily  implementable . The  proarams  could  be  rewritten 
without  these  limitations.  Whether  a revision  is 
justified  depends  on  two  factors:  1)  the  likelihood  of 

a similar  cruise  with  only  the  computing  power  of  the 
HP  9830  available,  and  2)  the  realization  that  the 
HP  9830  is  not  a general  purpose,  high  speed  computer. 
The  first  factor  will  certainly  be  met  many  times. 
Failure  to  appreciate  the  second  factor  could  lead  to 
unachievable  expectations. 


HP  9830  CAPABILITIES  AND  LIMITATIONS 


The  HP  9830  calculator  as  configured  on  the  LYNCH  had 
7904  words  of  memory,  a printer,  a 10"  x 15"  flat-bed 
plotter  and  one  cassette  read/write  unit.  The  calculator 
is  programmable  in  the  BASIC  language  with  many  additional 
features  to  allow  communication  with  the  cassette  and 
operating  system  of  the  calculator. 

The  calculator  keyboard  is  well  designed.  Keyboard 
entry  and  editing  is  easy  and  flexible.  Programming  in 
BASIC  is  simple  and  yet  quite  powerful.  There  is  a great 
deal  of  flexibility  in  the  use  of  the  cassette  for  proqram 
and  data  storage.  Indeed,  it  would  not  have  been  possible 
to  write,  debug  and  use  so  many  programs  in  such  a short 
time  if  the  HP  9830  were  not  such  an  easy  machine  to  use. 
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When  our  efforts  resulted  in  our  first,  neatly 
drawn  and  useful  chart,  we  were  very  pleased.  However, 
once  the  novelty  wore  off  and  the  task  of  data  entry, 
editing  and  charting  became  routine,  we  noticed  how 
slow  the  calculations  and  plotting  actually  were, 
compared  to  a general  purpose  minicomputer  and  graphic 
devices  The  calculator  does  not  have  large  amounts 
of  semiconductor  or  mass  storage  memory.  BASIC  programs 
are  not  precompiled  but  rather  recompiled  each  time  they 
are  run.  The  plotter  provided  an  effective  9 inch  by 
13.5  inch  plotting  area  inasmuch  as  a margin  had  to  be 
left  for  annotations.  We  were  able  to  draw  a grid 
30  minutes  in  latitude  and  50  minutes  in  longitude  at 
the  latitude  of  our  survey.  This  restriction  forced  us 
to  make  several  passes  through  the  data  before  we  had  a 
complete  set  of  charts. 


CONCLUSION 


An  underway  data  logging,  processing  and  graphic 
display  software  package  was  written  and  made  operational 
on  an  HP  9830  calculator  in  a short  time.  These  programs 
permitted  data  entry  and  merging  of  navigation  and 
magnetics  information,  editing  and  plotting  capabilities. 
Suitably  revised,  these  programs  could  be  used  on  other 
geophysical  cruises  where  only  limited  computing  power 
is  available.  Clearly  this  software  and  the  HP  9830 
calculator  cannot  compete  in  speed  or  flexibility  with  a 
larger,  general  purpose  minicomputer.  But  as  a basic, 
portable,  low  maintainance  hardware  and  software  system 
it  is  a very  attractive  package. 


6 


APPENDIX  I 


Hewlett  Packard  9830  Timing 


A simple  program  was  written  to  test  how  long  certain 
computational  operations  took.  The  FOR... NEXT  loop  was 
used  with  various  types  and  numbers  of  statements  within 
the  loop. 


# ITERATIONS 

WITHIN  LOOP 

TIME  (SECONDS) 

1 . 

1000 

no  instruction 

4 . 2 

2. 

1000 

10  REM 

12 

3. 

1000 

1 K1=K1+1 

12.5 

4 . 

1000 

10  K1=K1 + 1 

82 

5. 

1000 

1 Kl=Kl *1 

14 

6 . 

5000 

1 K1=K1 +1 

60.5 

7. 

10000 

1 Kl=Kl+ 1 

122 

Subtracting  line  3 from  line  4 above,  yields: 

1000*9  K1=K1+1  taking  69.5  seconds 
Or,  0.008  seconds  per  add  operation. 

Because  the  programs  are  executed  line  by  line,  you 
pay  a penalty  for  having  non  executing  statements  (like 
the  REMARK)  within  your  programs. 
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APPENDIX  2 


ALLOCATION  OF  COMMON  IN  PROGRAM 


The  variables  in  common  have  the  following  meaning: 

L - Array:  stores  information  about  available  space 

for  data  and  current  number  of  data  points  defined. 

L ( 1 ) - maximum  number  of  fixes  that  can  be 
processed 

L(2)  - maximum  number  of  data  points,  at 

even  five  minutes,  that  can  be  processed 
L(3)  - number  of  fixes  actually  in  memory.  The 
number  of  fixes  logged  plus  those  added 
by  the  edit  routine  minus  those  deleted 
L(4)  - number  of  even  five  minute  navigation 
points  available  in  memory 
L(5)  - number  of  defined  magnetic  anomaly  points, 
(less  than  or  equal  to  L(4)) 


Y - 

Array : 

string  array,  not  used. 

D - 

Array : 

the  day  number  array,  corresponding  to 
each  of  the  possible  L(3)  fixes  in  memory. 

T - 

Array : 

the  fix  time  array,  hour,  then  minute 
stored  separately. 

N - 

Array : 

the  fix  latitude  array,  degrees,  minutes. 

W - 

Array : 

the  fix  longitude  array,  degrees,  minutes. 

P - 

Array : 

the  five  minute  position  array,  latitude 
and  longitude,  in  degrees  and  fraction. 

A - 

Array : 

the  time  and  data  array. 

First  element:  even  five  minute  interpolated 

Second  element:  magnetic  anomaly  value. 

Third  element:  not  used. 
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APPENDIX  3 

PROGRAM  LOADING  FROM  CASSETTE 


With  the  correct  program  cassette  loaded  into  memory 
it  will  be  possible  to: 

1.  Log  navigation  data  into  memory  and  produce 
a listing  of  fixes,  speeds,  headings  and 
distances . 

2.  Edit  (insert,  delete,  replace  and  list) 
navigation  data. 

3.  Produce  five  minute  interpolated  positions 
and  times  from  the  navigation  data. 

4 . Calculate  the  total  magnetic  regional  field  and 
the  magnetic  anomaly  for  each  measured  magnetic 
field  value  entered  at  corresponding  five 
minute  interpolated  positions. 

5.  Edit  (delete,  replace  and  list)  magnetic  anomaly 
values . 

6.  Plot  a Mercator  chart  annotated  with  time. 

7.  Plot  a Mercator  chart  with  the  magnetic  anomaly 
plotted  in  profile,  parallel  to  the  ship's  track. 

The  directions  below  apply  to  cassette  #3.  This  cassette 
contains  programs  1-4  from  above.  These  programs  are  stored 
as  one  file.  REM  (comment)  statements  have  been  removed 
to  conserve  calculator  memory. 

LOADING  CASSETTE  #3: 

1.  Place  program  cassette  #3  into  the  cassette  reader 
and  close  the  door. 

2.  Type  in 

LOAD  KEY  0 


3.  Push  the  EXECUTE  key.  This  procedure  places  the 
correct  instructions  into  the  function  keys  fl-f4. 
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APPENDIX  3 (continued) 


4.  The  cassette  unit  will  start  reading  and  the 
display  will  blank  out.  When  the  cassette  has 

stopped  reading  the  'I ' symbol  will  again  appear  in 

the  display. 

5 . Type  in 

LOAD  1 

6.  Push  the  EXECUTE  key. 

7.  When  the  symbol  again  appears  in  the  display, 
indicating  that  the  cassette  unit  has  finished 
reading,  rewind  the  program  cassette  by  pushing 
the  REWIND  key.  This  key  is  located  in  the 
upper  right  hand  corner  of  the  keyboard. 

8.  After  the  cassette  unit  has  stopped  moving  (look 
into  the  cassette  unit  through  the  window), 
remove  the  cassette  from  the  read  unit. 

Four  programs  are  now  available.  They  are: 

1.  FIX  DATA  LOGGING  (fl) 

2.  FIX  EDIT  ( f 2 ) 

3.  FIX  INTERPOLATION  ( f 3 ) 

4.  MAGNETICS  DATA  LOGGING  AND  ANOMALY  CALCULATION  ( f 4 ) 

To  select  one  of  these  programs  push  the  correct 
function  key  located  on  the  upper  left  hand  part  of  the 
keyboard.  The  function  keys  are  given  in  parentheses  above. 
For  example,  to  run  the  fix  interpolation  program  push 
the  function  key  labeled  f3. 
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APPENDIX  4 


NAVIGATION  LOGGING  AND  FIX  CHECKING 


If  the  program  is  not  in  memory,  then  see  the 
detailed  instructions  (Appendix  3)  on  how  to  load  programs  from 
the  program  cassette.  Rewind  the  program  cassette 
by  pushing  the  'REWIND'  key  on  the  upper  right  corner 
of  the  keyboard.  Remove  the  cassette  after  it  has  been 
rewound . 

This  program  provides  the  means  to  enter  fixes  onto 
cassette  and  to  calculate  the  speed,  heading  and 
distance  between  fixes.  In  addition,  the  speed,  heading 
and  distance  calculations  can  be  done  for  data  already 
stored  on  cassette  or  already  stored  in  memory. 

Press  the  'fl'  function  key,  located  on  the  upper 
left  hand  corner  of  the  keyboard  to  begin  the  program. 

The  program  will  display  the  following  question: 

ENTER  FILE  # FOR  DATA  RESTORE? 

If  the  data  already  exists  on  cassette  and  you  want  a 
new  listing  of  speeds  and  headings,  merely  enter  the  file 
number  in  which  the  data  is  stored.  The  program  will  read 
the  data  from  the  cassette  into  memory  and  produce  a list 
of  the  fixes  with  the  calculated  speed,  heading  and 
distance  between  each  position.  The  program  will  go  to 
the  end  of  job  section. 

If  the  data  is  already  in  memory  or  if  you  want  to 
enter  a new  set  of  navigation  points,  type  in  a -1  in 
response  to  this  question. 


If  you  have  typed  in  a -1  the  program  then  asks: 
NEW  FIX  DATA? 


If  you  answer  "NO",  then  the  program  will 
already  in  memory  to  generate  the  list  of 
After  the  list  the  program  will  go  to  the 


use  the  positions 
speed  and  headings, 
end  of  job  section. 
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APPENDIX  4 (continued) 


To  enter  new  fix  information  respond  with  a "YES"  to 
the  question.  You  then  enter  the  day,  hour,  minute, 
latitude  degrees,  latitude  minutes,  longitude  degrees, 
longitude  minutes,  followed  by  the  END  OF  LINE  or 
EXECUTE  key.  North  latitude  and  west  longitude  are 
assumed  by  the  program  so  that  signs  (+  or  -)  need 
not  be  included.  For  example,  for  a fix  on  the  7th  at 
1026Z,  at  25°15.4'N,  72°37.8S'W  enter  the  following: 

7,  10,  26,  25,  15.4,  72,  37.85  (END  OF  LINE) 

After  the  last  fix  has  been  entered,  enter 

0,  0,  0,  0,  0,  0,  0 (END  OF  LINE) 

The  zero  entry  tells  the  program  that  there  are  no  more  fixes. 
A listing  of  the  fixes  with  speed,  heading  and  distance 
between  each  fix  will  be  printed.  After  the  listing 
has  been  produced,  the  program  goes  to  the  end  of  job 
section . 

End  of  job  section:  In  this  section  the  program  asks 

the  following  question: 

ENTER  FILE  # FOR  DATA  SAVING: 

If  you  would  like  to  save  these  fixes  on  cassette,  enter 
in  the  file  number  to  which  you  would  like  the  data  to  be 
written.  If  you  do  not  want  to  save  the  data  at  this 
step  enter  a -1.  Be  sure  the  data  cassette  is  ready  to 
accept  data  before  you  enter  the  file  number. 
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APPENDIX  5 


FIVE  MINUTE  INTERPOLATED  POSITIONS  BETWEEN  NAVIGATION  FIXES 


To  start  the  program  for  interpolating  positions 
between  fixes,  do  the  following: 

1.  Type  RUN  4000 

2.  Push  EXECUTE  key 

or,  hit  the  F3  function  key  if  the  calculator  has  been 
"set  up"  (Appendix  3)  for  this  option. 

Turn  on  the  printer  in  order  to  have  a copy  of  your 
inputs  and  for  the  output. 

The  program  has  two  options: 

1.  Returns  an  interpolated  latitude  and  longitude 
(degrees  and  minutes)  for  any  intermediate  day,  hour,  and 
minute  entered. 

2.  Lists  interpolated  positions  at  even  five  minute 
increments  between  fixes. 

The  second  option  is  the  one  used  for  most  navigation 
and  data  merging  in  this  suite  of  programs. 

Both  options  require  entering  the  two  positions  in 
response  to  successive  keyboard  quiries. 

To  enter  4 September  0641  25°33.2'N  069°16.1'W, 

type  4,  6,  41,  25,  33.2,  69,  16.1  (END  OF  LINE) 

Similarly  for  the  second  quiry,  to  enter 

4 September  1013  25°  38.1'N  068°57.4’W, 

type  4,  10,  13,  25,  38.1,  68,  57.4  (END  OF  LINE) 

The  quiry  "E  DAY  HR  MIN  OF  OBSERVATION"  will  appear  on 
the  console.  At  this  point  three  entries  can  be  made. 

1.  For  a specific  time  within  the  interval  between 
the  fixes,  other  than  the  time  of  the  first  fix,  enter  the 
day,  hour  and  time. 


> J 
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APPENDIX  5 (continued) 


For  example,  to  enter  4 September  0919, 

type  4,  9,  19  (END  OF  LINE) 

The  program  will  respond  with  the  interpolated 
position  and  pause.  To  make  another  entry,  hit  the 
space  bar  and  the  quiry  "E  DAY  HR  MIN  OF  OBSERVATION" 
will  be  repeated. 

2.  For  a listing  of  interpolated  positions  for  each 
even  five  minute  interval  between  fixes,  enter  the  day, 
hour  and  minute  of  the  first  fix.  In  the  example  above, 
type  4,  6,  41.  The  program  will  respond  with  a listing 
and  then  repeat  the  quiry  "E  DAY  HR  MIN  OF  OBSERVATION". 

3.  At  this  point,  if  you  wish  to  go  on  to  the  next 

navigation  fix,  enter  0,  0,  0 and  the  program  will  then 
move  the  second  fix  to  the  first  and  ask  for  a new 

second  fix.  If  you  wish,  instead,  to  enter  two  other  fixes 
type  STOP,  EXECUTE,  RUN4000,  EXECUTE  to  restart  the  program 


NOTE: 


The  preceeding  directions  do  not  apply  to  the  inter- 
polating program  used  with  the  overall  processing 
software  package.  This  package  allows  for  data  to  be 
saved  in  memory  and  on  cassette  for  later  use  by  other 
programs . 

To  use  this  program  within  the  processing  package, 
proceed  as  follows: 

1.  Start  the  program. 

2.  The  program  will  request  the  cassette  file  number 
from  which  navigation  data  can  be  read  back  into 
memory.  Type  in  the  file  number  for  the  data 
file  desired.  If  the  data  is  already  in  memory, 
then  enter  a -1  (minus  one)  . 

3.  The  program  will  calculate  even  five  minute 
interpolated  times  and  positions  between  each  fix. 
A listing  of  these  will  be  made  at  the  same  time. 
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APPENDIX  5 (continued) 


4.  After  all  calculations,  the  program  will 
request  the  file  number  on  which  you  would 
like  the  data  stored.  Usually,  this  number 
will  be  the  same  as  the  input  file  number  entered 
at  step  2.  If  you  do  not  want  to  save  the  data 
at  this  time,  enter  a -1  (minus  one). 
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APPENDIX  6 


NAVIGATION  EDIT 


This  routine  allows  you  to  edit  the  day,  hour, 
minute,  latitude  and  longitude  of  every  fix  stored 
in  memory.  New  fixes  can  be  inserted  and  old  fixes 
can  be  deleted.  The  basis  for  editing  any  fix  is  its 
key  number.  This  number  corresponds  to  the  fix's  record 
number  (or  position)  in  storage.  Since  there  can  be  up  to 
thirty  fixes  in  memory  at  one  time,  there  are  up  to  thirty 
keys.  Key  7,  for  example,  corresponds  to  fix  7.  To  help 
determine  the  key  number  of  a particular  fix,  two  commands 
are  available.  The  L,  or  list,  command  produces  a list 
of  selected  fixes  from  memory,  along  with  their  keys. 

The  F,  or  find,  command  locates  a particular  fix  and  key 
number  by  day,  hour  and  minute. 

Initially  the  program  requests  that  you  enter  the 
cassette  file  number  containing  the  data  to  be  edited. 

The  prompt  is: 

ENTER  FILE  # FOR  DATA  RESTORE? 

Type  in  the  cassette  file  number  where  your  data  is 
located.  If  the  data  is  already  in  memory,  type  in  a 
-1  (minus  one).  If  you  have  specified  a valid  file 
number,  the  cassette  will  be  positioned  and  the  data  read 
into  memory. 

The  program  will  next  print: 

# OF  FIXES  IN  ARRAY  = NN 

# OF  INTERPOLATED  TIME  POINTS  IN  ARRAY  = MM 

where  values  of  NN  and  MM  will  depend  on  the  data  file 
chosen.  NN  should  be  between  1 and  30.  MM  should  be 
between  0 and  200. 

The  program  prompts  for  edit  commands  with  the  line: 

EDIT  OPTION  (L, F, R, D, I , E) ? 


1 


r 
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APPENDIX  6 (continued) 


This  line  is  referred  to  as  the  'major  options  list'. 
The  meaning  and  further  input  for  each  option  is 
described  below. 


LIST  (L) 


The  list  option  provides  a means  to  list  part  or  all 
of  the  navigation  data,  in  sequential  order,  with 
corresponding  key  numbers  for  each  fix.  The  LIST  option 
prompts  with: 

LIST  FROM-TO  KEYS? 

Type  in  the  first  and  last,  inclusive,  keys  that  you  want 
listed.  If  the  first  number  specified  is  less  than  zero, 
then  you  will  return  to  the  major  options  list.  If  the 
second  key  number  is  larger  than  the  number  of  fixes  in 
memory,  it  will  be  redefined  as  the  last  fix.  The  first  key 
number  must  be  less  than  or  equal  to  the  second  number. 


FIND  (F) 


The  find  option  allows  you  to  determine  the  key  number 
of  any  specific  fix  in  memory.  The  find  option  prompts 
with: 

FIND:  DAY,  [(OUR, MINUTE? 

Enter  the  day,  hour  and  minute  of  the  fix  you  are  looking 
for.  The  program  will  print  the  key  number  for  this  fix  if 
it  is  located.  If  it  cannot  find  the  fix  it  will  display 
the  message: 

DATE  NOT  FOUND  day  time 

where  the  'day'  and  'time'  are  the  day  and  time  searched 
for . 

After  printing  the  key  number  or  displaying  the  above 
message  you  return  to  the  major  options  list.  If  you  enter 
a negative  value  for  the  day,  the  program  returns  immediately 
to  the  major  options  list. 
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APPENDIX  6 (continued) 


REPLACE  (R) 


The  replace  option  allows  you  to  modify  a fix 
already  stored  in  memory.  The  replace  option  asks 
for  the  key  number  of  the  fix  to  be  replaced: 

REPLACE  KEY  #? 

Enter  the  key  number  of  the  fix  you  want  to  modify. 
Entering  a zero  or  negative  number  will  cause  the 
program  to  return  to  the  major  options  list. 

Once  you  have  selected  a key  number,  the  program 
prompts  you  by  displaying,  in  turn,  the  current  value 
for  day,  hour,  minute,  latitude  and  longitude.  As  each 
current  (or  old)  value  is  displayed,  you  enter  in  the  new 
value.  The  program  assumes  west  longitude  so  the  sign  of 
the  longitude  need  not  be  entered.  It  will  be  made 
negative  automatically.  No  check  is  made  for  valid  data 
or  backwards  in  times.  Use  the  fix  logging  routine  after 
editing  any  navigation  data  to  check  for  correct  positions 
and  times. 


DELETE  (D) 


The  delete  option  allows  you  to  delete  an  existing 
fix  from  memory.  It  prompts  with: 

DELETE  KEY  #? 

Enter  the  key  number  to  delete.  A number  less  than  or 
equal  to  zero  will  return  you  to  the  major  options  list. 
After  the  record  is  deleted  the  program  prints: 

KEY  # n DELETED 

and  returns  to  the  major  options  list.  If  *n'  is  greater 
than  the  number  of  fixes  in  memory  you  will  be  asked  to 
enter  the  key  number  again. 

All  records  after  the  record  deleted  are  renumbered. 
That  is,  if  key  number  3 is  deleted,  the  record  associated 
with  key  number  4 now  is  now  number  3,  etc. 
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APPENDIX  6 (continued) 


INSERT  ( I ) 


The  insert  option  allows  you  to  add  new  navigation 
points  to  the  data  already  stored  in  memory.  The 
insert  option  prompts  with: 

INSERT  DATA  BEFORE  KEY  #? 

Type  in  the  key  number  of  the  record  before  which  you 
want  the  new  fix  to  be  placed.  Entering  a 0 or  negative 
number  returns  you  to  the  major  options  list. 

After  typing  in  a valid  key  number  the  program 
prompts  with: 

INSERT  DAY,  HR,  MIN,  LAT,  LONG? 

Enter  in  the  new  fix  day,  hour,  minute,  latitude  (degrees, 
minutes)  and  longitude  (degrees,  minutes).  West  longitude 
is  assumed.  This  is  the  same  form  as  used  by  the  fix 
logging  program.  You  need  not  enter  the  minus  signs  for 
west  longitude.  No  check  is  made  for  valid  positions, 
or  times.  Run  the  fix  logging  program  to  check  speeds 
and  headings  after  any  edit  to  the  navigation  data. 

All  records  after  the  newly  added  fix  are  renumbered 
with  new  keys.  The  new  key  is  the  old  key  plus  one. 


END  (E) 


When  you  are  done  editing,  type  in  the  "E"  command 
to  go  to  the  end  of  job.  The  program  will  prompt  with: 

SAVE  DATA  BEFORE  LOGGING!!! 

ENTER  FILE  # FOR  DATA  SAVING? 

Type  in  the  cassette  file  number  on  which  you  want  to 
save  the  edited  data.  If  you  do  not  want  to  save  the 
data  at  this  time,  enter  in  a -1  (minus  one) . 


FIX  EDIT  ROUTINE 
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APPENDIX  7 


OBSERVED  MAGNETIC  FIELD  LOGGING  AND  ANOMALY  CALCULATION 


If  the  programs  are  not  in  memory,  then  see  the 
detailed  instructions  on  how  to  load  programs  from 
a program  cassette  (Appendix  3) . 

A bi-directional  linear  interpolation  scheme  is  used 
to  calculate  the  regional  field  value  within  the  degree 
square  of  interest.  The  even  five  minute  interpolation 
program  must  be  run  prior  to  this  program. 

To  start  the  program  press  the  * f4'  function  key 
located  on  the  upper  left  hand  corner  of  the  keyboard. 

The  program  will  display  the  following  question: 

ENTER  FILE  # FOR  DATA  RESTORE? 

To  load  the  navigation  data  from  cassette,  have  the 
cassette  in  the  cassette  unit,  and  enter  the  file  number 
for  the  correct  data  as: 

4 (END  OF  LINE) 

This  will  cause  the  data  from  file  4 to  be  read  into 
memory.  If  the  navigation  data  is  already  in  memory 
then  enter  a -1  (minus  one)  as: 

-1  (END  OF  LINE) 

Next  the  program  requests  that  you  enter  the 
START  POINTER.  This  option  allows  you  to  start  processing 
at  other  than  the  first  time  and  position  of  the  data. 

To  start  at  the  first  data  point  enter  a 1 as : 

1 (END  OF  LINE) 

To  start  at  the  57-th  point  enter  a 57  as: 

57  (END  OF  LINE) 

When  the  program  needs  new  values  for  the  regional 
field  it  will  print  the  current  latitude  and  longitude  in 
the  new  degree  square.  Then  it  requests  that  you  enter 
the  degree  bounds  which  surround  this  data  point  and  the 
corresponding  regional  field  values  for  the  corners  of 
this  degree  square.  For  example,  the  program  miqht 
display  the  following  line. 
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APPENDIX  7 (continued) 


DATA  AT  25.001  -70.537  (note:  West  is  negative) 

ENTER  THE  DEGREE  BOUNDS 

The  degree  bounds  which  surround  this  data  point  are: 

Top:  26°N;  Bottom  25°N;  Left:  71°W;  and  Right:  70°W. 

The  numbers  are  entered  in  this  order: 

26,  25,  71,  70 

Since  north  latitude  and  west  longitude  are  assumed,  no 
signs  (+  or  -)  need  to  be  included.  The  order  is  top, 
bottom,  left,  right. 

The  regional  field  values  for  each  square  for  the  time 
of  the  cruise  has  been  calculated  and  tabulated.  The  table 
gives  the  following  information: 

DEGREE  CROSSING  FIELD  VALUE  (gammas) 

26  N,  71  W 48137 

25  N,  71  W 47497 

25  N,  70  W 47338 

26  N,  70  W 47973 

The  computer  will  request  the  field  values  by  typing: 

ENTER  4 REG.  FIELD  VALUES 

Enter  the  values  in  the  order  given  above.  That  is,  top 
left,  bottom  left,  bottom  right,  top  right: 

48137,  47497,  47338,  47973  (END  OF  LINE) 

For  every  data  point  within  the  degree  square,  the 
program  will  display  the  time  as 

ENTER  VALUE  FOR  1835? 

Type  in  the  measured  value  of  the  total  magnetic  field 
for  this  time.  If  there  is  no  value,  then  enter  a 0 
(zero).  For  example,  if  the  measured  value  is  47332  gammas, 
enter 


47332  (END  OF  LINE) 
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APPENDIX  7 (continued) 


If  you  do  not  want  to  enter  any  more  values,  then  enter 
-1  (minus  one) . Any  times  skipped  will  have  the 
anomaly  value  set  to  30000.  Any  measured  value  larger 
than  60000  gammas  or  smaller  than  30000  gammas  will  not 
be  accepted  and  the  program  will  ask  for  the  value  again 

When  there  is  no  more  data  or  if  you  have  given  a 
-1  as  the  value  for  the  measured  field,  the  program  will 
ask : 

DO  YOU  WANT  A LIST  OF  THE  DATA? 

If  you  do,  then  respond  with  a "YES".  If  you  do  not 
want  a list,  then  answer  "NO". 

Finally,  the  program  asks: 

ENTER  FILE  # FOR  DATA  STORAGE: 

To  save  the  results  on  cassette,  enter  the  file  number 
of  the  cassette  at  which  you  want  the  data  saved.  For 
example,  to  save  the  data  on  file  7 enter 

7 (END  OF  LINE) 

Ordinarily  the  data  is  read  back  onto  the  same  cassette 
file  as  it  was  read  from.  If  you  do  not  want  to  store 
the  data  on  cassette,  then  enter  a -1  (minus  one)  as 


-1  (END  OF  LINE) 


MAGNETIC  LOGGING  AND  ANOMALY  CALCULATION 
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APPENDIX  8 


MAGNETICS  EDIT  ROUTINE 


The  intent  of  the  magnetics  editing  routine  is  to 
locate  and  correct  erroneous  geomagnetic  data  Magnetic 
data  will  have  already  been  logged  and  magnetic  anomaly 
values  calculated  via  program  MAGNETIC  LOGGING  AND  ANOMALY 
CALCULATION  (Appendix  7) . Three  methods  are  used  to  detect 
bad  values.  These  are  1)  a smoothness  criterian  (excessive 
second  drivative) , 2)  unreasonable  value  checking  (lower 
limit),  and  3)  excessive  deviation  from  median  value  of 
three  points. 

If  the  program  is  not  in  memory,  than  see  the  detailed 
instructions  (Appendix  3)  on  how  to  load  programs  from  the 
program  cassette.  Rewind  the  program  cassette  by  pushing 
the  'REWIND'  key  on  the  upper  right  corner  of  the  keyboard. 
Remove  the  cassette  after  it  has  been  rewound.  The 
printer  must  be  on  for  this  program. 

The  program  will  display  the  following: 

MAGNETIC  EDTT  ROUTINE 

ENTER  FILE  ft  FOR  DATA  RESTORE 

Enter  the  file  number  of  the  magnetics  data  that  you  want 
placed  into  memory.  The  data  cassette  should  be  in  the 
cassette  reader.  If  the  data  is  already  in  memory  then 
enter  -1.  If  you  have  indicated  a file  number  greater  than 
or  equal  to  zero,  the  data  will  be  road  from  the  cassette. 
The  program  will  print: 


FILE  ATTRIBUTES 


L ( 3 ) = X L ( 4 ) = Y 


L ( 5)  = 7 


where  X is  the  number  of  valid  fixes  in  memory,  Y is  the 
maximum  number  of  even  five  minute  navigation  (and  data) 
points  available  and  7.  is  the  number  of  defined  magnetic 
anomaly  points. 


The  program  will  print 


NEW  VALUE  FOR  L ( 5 ) 
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Enter  the  desired  number  of  defined  magnetic  anomaly  points. 
Entering  a value  of  zero  or  less  will  leave  L(5)  unchanged. 
The  purpose  of  this  option  is  to  permit  deleting  of 
erroneous  (undefined)  values  of  magnetic  anomaly  points 
at  the  end  of  the  data  array.  The  magnetics  logging  program 
can  then  be  rerun  starting  with  this  new  value  of  L(5). 

The  program  will  print 

NOW  L ( 5 ) = W 

to  show  that  the  value  of  L(5)  is  now  equal  to  W. 


Next  the  program  executes  the  standard  purge  section 
where  all  anomaly  values  less  than  or  equal  to  ten  gammas 
are  redefined  to  30000,  the  null  or  'no  data'  indication. 

Ten  gammas  is  a reasonable  cut-off  value  for  this  survey 
area.  The  program  prints 

n VALUES  REDEFINED 

to  show  how  many  points  exceeded  the  cut-off  value. 

The  program  then  executes  the  smoothness  check  section. 
The  message 

SMOOTHNESS  CHECK  SECTION 

2-nd  DERIVATIVE  MAXIMUM  = n % OF 
MEDIAN  ALLOWED  = m 

is  printed  and  each  magnetic  anomaly  point  in  turn  is 
checked  for  smoothness  (second  derivative)  and  excessive 
deviation  from  median  of  surrounding  points.  The  values 
of  n and  m are  set  to  35  gammas/kilometer  and  10  percent, 
respectively.  When  a value  does  not  pass  the  smoothness 
check  the  following  is  printed: 

SMOOTHNESS  ERROR 

KEY  = n PAST  = X PRESENT  = Y FUTURE  = Z 


TIME 


T 


2-nd  = W 


•»  ■ 
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where  n refers  to  the  key  number  of  the  present  value; 

X,  Y,  and  Z refer  to  the  past,  present  and  future  values 
of  the  magnetic  anomaly  data;  T is  the  time  of  the  present 
(or  current  value  under  consideration)  and  W is  the  cal- 
culated value  of  the  second  derivative.  When  a value  does 
not  pass  the  smoothness  check  it  will  then  be  checked  for 
excessive  deviation  from  median  value.  If  this  check  is 
not  passed  then  the  following  is  printed; 

BAD  MEDIAN,  VALUE  = Y KEY  = n 

where  Y is  the  anomaly  value  and  n is  its  key  number. 

Many  values  which  only  failed  the  smoothness  check  were 
found  to  be,  in  fact,  valid  values.  However,  most  values 
which  failed  both  tests  were  found  to  be  incorrect  on 
further  examination.  After  all  points  in  the  array  have 
been  checked,  the  following  message  will  be  printed 

p SMOOTHNESS  ERRORS  FOUND 

where  p is  the  number  of  points  failing  either  or  both  tests. 
The  program  will  display 

MAGNETIC  KEY  # TO  REDEFINE? 

Enter  in  the  key  number  for  the  magnetic  value  you  wish 
to  redefine.  Consult  the  results  of  the  smoothness  check 
section  in  order  to  select  the  key  number.  The  program 
will  print 

OLD  VALUE  = Y NEW  = ? 

Enter  in  the  new  value  you  would  like  to  assign.  Y is 
the  original  value  in  gammas.  The  program  will  continue 
requesting  a key  number  unt i 1 a zero  or  negative  value  is 
supplied  for  the  key  number.  The  program  then  requests 

START  - END  KEYS  FOR  LIST? 

Enter  the  first  and  last  key  number  of  records  you 
would  like  to  list  out  on  the  printer.  These  records 
are  printed  under  the  following  heading: 

#,  HOUR,  LATITUDE,  LONGITUDE,  and  ANOMALY 
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The  start  day  (START  DAY  = ) is  printed  at  the  start  of  the 
list  for  reference.  Enter  negative  values  for  no  list. 

Finally,  the  program  gives  you  the  option  of  storing 
the  modified  data  back  onto  cassette.  The  program  will 
display : 

ENTER  FILE  # FOR  DATA  STORAGE? 

Enter  a negative  number  if  you  do  not  want  to  write  the 
data  back  to  cassette.  The  program  ends  with 


END  OF  MAGNETICS  EDIT. 


I1AGNETICS  EDIT  ROUTINE 
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APPENDIX  9 


NAVIGATION  PLOTTING  PROGRAM 


This  program  plots  navigational  fixes  on  an  operator 
specified  Mercator  grid  and  annotates  the  time  at  right 
angles  to  the  track  line. 

To  load  the  program,  insert  the  program  cassette, 
and  load  the  file  by  typing  LOAD  < file  # >.  Start 
the  program  by  typing  "RUN".  "REWIND"  the  program 
cassette  and  insert  the  data  cassette. 

The  operator  is  queried  five  times  as  follows: 

1.  Mercator  grid  parameters.  Enter  the  longitude 
scale  in  inches  per  degree,  and  the  number  of  minutes 
between  grid  lines  desired.  Thirty  minutes  would  draw 
two  lines  for  each  degree. 

2.  File  number  of  data  to  be  plotted  (enter  -1,  if 
data  are  already  in  core) ; option  to  list  the  positions 
plotted  on  the  printer  (1  to  list,  -1  not  to);  minimum 
spacing  between  annotations  in  inches  (prevents 

successive  number  of  annotations  from  overprinting) . The  minimum 
separation  should  be  0.1  inch,  the  size  of  the 
numbers.  On  a smaller  scale  chart,  a larger  value  is 
suitable.  After  this  entry,  there  will  be  a pause  while 
the  data  file  is  loaded  into  core. 

3.  Bounds  of  the  plot  in  degrees  and  minutes.  Enter 
in  the  following  order  - top,  bottom,  left,  right. 

4.  The  terminal  will  display  "SET  PLOTTER"  and  then 
give  the  dimensions  to  set  the  plotter  in  inches.  These 
dimensions  include  a one  inch  border  on  the  left,  and  a half  inch 
border  on  the  top,  bottom  and  right  on  the  10"xl5"  plotter.  After 
setting  the  plotter  pen  to  the  specified  dimensions  using 

the  LOWER  LEFT  and  UPPER  RIGHT  controls,  resume  the  proqram 
by  typing  "CONT"  and  "EXECUTE". 
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APPENDIX  9 (continued) 


5.  The  final  query  is  a label  for  the  left-hand 
margin  of  the  plot.  Make  the  entry  slowly  because  the 
plotter  annotates  each  letter  as  you  type  it.  Wait  for 
"ENTER  LABEL  AND  RUN  TIME"  to  reappear  before  typing 
the  next  letter  or  space.  When  finished  with  the  last 
letter,  type  the  "STOP"  key  and  the  program  will  resume. 
When  circumstances  might  involve  making  repeated  plots 
of  the  same  fixes,  such  as  in  editing  the  navigation, 
it  is  useful  to  include  the  run  time  and  date  in 
the  label,  so  to  facilitate  identification  of  the 
chronology  of  the  plots. 

The  grid  is  drawn  and  the  bounds  labelled.  The 
positions  are  marked  with  a dot  and  the  time  annotated. 

If  the  list  option  in  entry  2 is  taken,  a list  of  all  the 
fix  positions  and  times  falling  within  the  specified  grid 
will  be  printed.  At  the  conclusion  of  the  plotting,  the 
printer  will  log  the  number  of  points  and  display  "END  OF 
DATA"  on  the  calculator. 


To  plot  additional  files  on  the  same  grid,  type  "CONT" 
and  "EXECUTE".  The  terminal  will  ask  for  the  number  of  the 
next  file  on  the  data  cassette  (which  could  be  on  another 
cassette) . Be  sure  to  "REWIND"  each  cassette  before 
removing  it  from  the  reader,  to  reduce  the  chance  of 
tape  damage. 
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APPENDIX  10 


PROGRAM  FOR  PLOTTING  DATA  ALONG  SHIP'S  TRACK 


This  program  plots  data  along  the  ship's  track  using 
the  track  line  as  a datum.  The  perpendicular  displacement 
from  the  track  is  proportional  to  an  operator  specified 
scale.  The  track  is  plotted  on  a Mercator  grid  as  a 
series  of  dots  and  the  data  as  a continuous  trace.  The 
trace  is  broken  at  course  changes  of  more  than  thirty 
degrees  and  if  there  is  a gap  in  the  data.  The  track 
line  can  be  assigned  an  arbitrary  value  so  that  the 
fluctuations  in  the  plotted  parameter  can  be  displayed 
at  a larger  scale  if  desired.  Positive  values  are  plotted 
on  the  "top"  side  of  the  track,  minus  values  below. 

To  load  the  program,  insert  the  program  cassette 
and  load  the  file  by  typing  LOAD  < file  # >.  Start  the 
program  by  typing  "RUN".  "REWIND"  the  program  cassette 
and  insert  the  data  cassette. 

The  operator  is  queried  six  times  as  follows: 

1.  Mercator  grid  parameters.  Enter  the  longitude 
scale  in  inches  per  degree,  and  the  number  of  minutes 
between  grid  lines  desired.  Fifteen  minutes  would  draw 
four  lines  for  each  degree. 

2.  File  number  of  data  cassette  to  be  plotted 
(enter  -1,  if  data  is  already  in  core);  option  to  list  the 
time,  position  and  data  plotted  (enter  "1"  for  listing, 

"-1"  not  to) . After  this  entry  there  will  be  a pause  while 
the  data  are  loaded  from  the  specified  cassette  file. 

3.  Two  parameters  controlling  the  deflection  and 
placement  of  the  data  track  in  respect  to  the  track  line 
are  requested.  First  the  scale  in  scientific  units  per 
inch  displacement  from  the  track;  second  a constant  bias 
value  to  be  subtracted  from  all  data  values.  For  plotting 
the  geomagnetic  field  anomaly,  500,  0 are  typically  entries. 
For  bathymetry  -200,  2000  might  be  entered,  where  2000  is 
the  average  depth  of  the  water.  The  minus  200  is  used 

to  invert  the  sounding  profile  so  that  shallower  values 
plotted  above  deeper  values,  since  soundings  are  all 
negative  in  respect  to  the  sea  surface. 
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APPENDIX  10 (continued) 


4.  Bounds  of  the  plot  in  deqrees  and  minutes. 

Enter  in  the  following  order  - top,  bottom,  left,  right. 

Latitude  north  is  positive,  and  longitude  west  is 
negative  in  the  program.  If  the  minus  sign  is  not 
entered  the  plot  will  be  reversed  right  for  left. 

5.  The  calculator  will  display  "SET  PLOTTER"  and  then 
give  the  dimensions  to  set  the  plotter  in  inches.  These 
dimensions  include  a ore  inch  border  on  the  left,  and 

a half  inch  border  on  the  top,  right  and  bottom.  After 
setting  the  plotter  pen  to  the  specified  dimensions 
using  the  LOWER  LEFT  and  UPPER  RIGHT  controls,  resume 
the  program  by  typing  "CONT"  and  "EXECUTE". 

6.  The  final  query  is  a label  for  the  left-hand 
margin  of  the  plot.  Make  the  entry  slowly  because  the 
plotter  annotates  each  letter  as  you  type  it.  Wait  for 
"ENTER  LABEL  AND  RUN  TIME"  to  reappear  before  typing  the 
next  letter  or  space.  The  label  should  include  the 

scale  and  bias,  as  well  as  the  file  description  and  date  of  run. 
The  entire  left  margin  can  be  used.  When  finished  type  the 
"STOP"  key  and  the  plotting  will  begin. 

The  grid  is  drawn  first  and  the  bounds  labelled.  The 
ship's  track  is  indicated  by  a series  of  dots  as  the  program 
scans  the  data  file.  The  data  points  are  then  plotted 
by  a smooth  trace  as  the  program  scans  the  data  file  a 
second  time.  If  the  list  opti  n in  entry  2 is  opted, 
a list  of  all  the  times,  positions  and  data  values  within 
the  grid  is  given  on  the  printer.  At  the  end  of  the  data 
file  the  program  will  log  the  number  of  data  points 
plotted . 

To  plot  additional  files  on  the  same  grid,  type  "CONT" 
and  "EXECUTE".  The  terminal  will  ask  for  the  number  of  the 
next  file  on  the  data  cassette  (which  could  be  on  another 
cassette) . Be  sure  to  "REWIND"  each  cassette  before 
removing  it  from  the  reader,  to  reduce  the  chance  of 
tape  damage. 
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